Vibrational Spectra of Small Protonated Peptides from Finite Temperature MD Simulations and IRMPD Spectroscopy.

Finite temperature Born-Oppenheimer DFT-based molecular dynamics simulations are presented for the vibrational spectroscopy of the prototype gas-phase Ala2H(+) and Ala3H(+) protonated peptides. The dynamics and the vibrational signatures are used to interpret IR-MPD spectra recorded in the NH/OH stretch region. Molecular dynamics simulations are one way to go beyond the harmonic approximations commonly applied for the calculations of infrared spectra, naturally including all anharmonicities, i.e. mode couplings, vibrational and dipole anharmonicities. The dynamics of the peptides allows understanding of the evolution of the shape and width of the N-H bands when increasing the size of the peptide, as demonstrated here with the two small prototypes Ala2H(+) and Ala3H(+). Hence, the conformational dynamics of Ala2(+) at room temperature participates to the broadening of the IR active bands. The complex N-H broadband of Ala3H(+) is shown to result from the dynamics of the N-H groups in the different peptide families, with a special role from breaking/reforming of hydrogen bonds involving the N-H groups. Taking this dynamics into account is thus mandatory for the understanding of this band in the 300-400 K experimental spectrum.

Molecular Dynamics and Room Temperature Vibrational Properties of Deprotonated Phosphorylated Serine.

The local structure of phosphorylated residues in peptides and proteins may have a decisive role on their functional properties. Recent IRMPD experiments have started to provide spectroscopic signatures of such structural details; however, a proper modeling of these signatures beyond the harmonic approximation, taking into account temperature and entropic effects, is still lacking. In order to bridge this gap, DFT-based Car-Parrinello molecular dynamics simulations have been carried out for the first time on a phosphorylated amino acid, gaseous deprotonated phosphoserine. It is found that all vibrational signatures are successfully reproduced, and new deconvolution techniques enable the assignment of the vibrational spectrum directly from the dynamics results and the comparison of vibrational modes at several temperatures. The lowest energy structure is found to involve a strong hydrogen bond between the deprotonated phosphate and the acid with relatively small free energy barriers to proton transfer; however, we find that proton shuttling between the two sites does not occur frequently. Anharmonicities turn out to be important to reproduce the frequencies and shapes of several experimental bands. Comparison of room temperature and 13 K, effectively harmonic dynamics, allows insight to be obtained into vibrational anharmonicities. In particular, a significant blue-shift and broadening of the C═O stretching frequency from 13 to 300 K can be ascribed to intrinsic anharmonicity rather than to anharmonic coupling to other modes. On the other hand, significant couplings are found for the stretching motions of the hydrogen bonded P-O bond and of the free P-OH bond, mainly with modes within the phosphate group.